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Abstract 

Bose-condensed gases are considered with an effective interaction strength varying 
in the whole range of the values between zero and infinity. The consideration is based on 
the usage of a representative statistical ensemble for Bose systems with broken global 
gauge symmetry. Practical calculations are illustrated for a uniform Bose gas at zero 
temperature, employing a self-consistent mean-field theory, which is both conserving 
and gapless. 
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1 Introduction 



The properties of systems with Bose-Einstein condensate are currently a topic of great inter- 
est, both experimentally and theoretically (see review articles [1-7]). One usually considers 
weakly interacting Bose gases, whose theory was pioneered by Bogolubov [8,9]. Binary 
atomic interactions in such gases can be modelled by contact potentials expressed through 
effective scattering lengths. But the latter can also be made rather large by means of 
the Feshbach resonance technique, so that effective atomic interactions could become quite 
strong [7,10,11]. Extension of the Bogolubov theory to Bose systems with strong interactions 
confronts the well known problem of conserving versus gapless approximations, as was for- 
mulated by Hohenberg and Martin [12]. This dilemma has recently been discussed in detail 
in the review paper by Andersen [5]. 

The Hohenberg-Martin dilemma of conserving versus gapless theories can be resolved by 
employing representative statistical ensembles [13]. Using such an ensemble for Bose systems 
with broken global gauge symmetry makes it straightforward to get a self-consistent theory, 
both conserving as well as gapless in any given approximation. In particular, the Hartree- 
Fock-Bogolubov (HFB) approximation, which is by construction conserving, can also be 
made gapless [14]. 

In the present paper, we consider an equilibrium Bose system with Bose-Einstein conden- 
sate. The main new results are twofold. First, we give a general mathematical foundation for 
the construction of the grand Hamiltonian for an arbitrary equilibrium system with broken 
gauge symmetry. The derivation of the grand Hamiltonian and the corresponding equations 
of motion are valid for any Bose system, whether uniform or nonuniform. Second, in the 
frame of a self-consistent mean-field theory for a uniform Bose gas, we study, both analyti- 
cally and numerically, the zero-temperature characteristics as functions of the gas parameter, 
varying the latter from zero to infinity. Specifically, the condensate fraction, sound velocity, 
normal and anomalous averages, and the ground-state energy as functions of the gas param- 
eter are investigated. The results are in good agreement with available computer Monter 
Carlo simulations. 

We use the system of units, where h = 1 and k B = \. 

2 Representative Ensemble for Bose-Condensed Sys- 
tems 

The description of a spinless Bose system at temperature T > T c above the condensation 
temperature T c can be done in terms of the field operators ip(r,t) and i{)'(r,t) dependeing 
on the spatial vector r and time t. The operators from the algebra of observables and 
other physical operators are defined in the Fock space J-{ij}) generated by the field operator 
ift . The related mathematical details of constructing the Fock space can be found in 

books [15,16]. Under the total number of particles N =< N >, being the average of the 
number-of-particle operator N, the grand Hamiltonian has the standard form 

H[^]=H[^]- l iN (T>T C ), 

where H[ip] is the Hamiltonian energy, which is invariant under the global gauge transfor- 
mations from the U(l) symmetry group. 
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At temperatures T < T c , the global gauge symmetry becomes broken. This is achieved 
by means of the Bogolubov shift [17,18] for the field operators 

^(r, t) — > ^ (r, t) = V (r, t) + ^i(r, t) , (1) 

in which r)(r,t) is the condensate wave function and ipi(r,t) is the field operator of uncon- 
densed atoms, enjoying the same Bose commutation relations as ip. The condensate wave 
function r)(r,t) is the system order parameter. Now all operators of physical quantities are 
defined on the Fock space generated by the field operator It is important to 

emphasize that the Fock spaces and are mutually orthogonal [13,19]. 

Thus, below T c , instead of one operator variable tp, there appear two variables, rj and ipi. 
These are linearly independent, being orthogonal to each other, 

J V *(r,t)Mr,t)dr = 0. (2) 

For two linearly independent variables, there are two normalization conditions. One is the 
normalization of the condensate function to the number of condensed atoms 

N = J \ V (r,t)\ 2 dr . (3) 

And another normalization condition is for the operator 

A>! = J 4(T,t)Mr,t) dr , (4) 

whose average yields the number of uncondensed atoms 

Ni = < A > . (5) 

The normalization condition (3) can be represented in the same form of the statistical average 
(5) by using the operator 

N = AT i^ , 

where 1^- is the unity operator in Fijpi). Then Eq. (3) is equivalent to the normalization 

A = < A > . (6) 

The statistical average of an operator A is defined in the standard way as 

< A > ee Tr p A , 

where p is a statistical operator and the trace is over Fijpi). 
One more restriction is 

<Vi(r,0>= 0, (7) 

which guarantees the conservation of quantum numbers. This can also be rewritten as the 
quantum conservation condition 

< A > = (8) 

for the self-adjoint operator 



A ee f [A(r, i)V>I(r, t) + A*(r, t)^(r, *)] dr , (9) 



in which X(r,t) is a complex function. 

Two other common conditions is the normalization of the statistical operator p, 

<ir>= 1 , (10) 

and the definition of the internal energy 

E =<H> (11) 

as the average of the Hamiltonian energy operator H = H[ip], which is a functional of the 
shifted field operator (1). 

An equilibrium statistical ensemble for a Bose-condensed system is the pair {^(ipi), p} 
of the space of microstates J 7 (ipi) and a statistical operator p. The notion of a representative 
ensemble stems from the works of Gibbs [20] , who emphasized that for the correct descrip- 
tion of the given statistical system, in addition to the standard conditions (10) and (11), it 
is necessary to take into account all other constraints that uniquely define the considered 
system. The corresponding statistical operator can be found from the maximization of the 
Gibbs entropy S = — Trplnp under the given statistical conditions. The conditional maxi- 
mization of the entropy is equivalent to the unconditional minimization of the information 
functional [16,21]. In the present case, in addition to conditions (10) and (11), we must also 
take into account the normalization conditions (5) and (6) and the conservation constraint 
(8). Hence the information functional is 

I[p]=Trp lnp + A (Trp- 1) + 

+(3(Tr pH-E)- (3p (Tr p A - N ) - (3^(Tr p N\ — N x ) — f3Ti p A , (12) 

in which Ao, f3, (3po, Pl^i, and (5\ are the appropriate Lagrange multipliers guaranteeng the 
validity of conditions (5), (6), (8), (10), and (11). Minimizing functional (12), we get the 
statistical operator 

P = \ e-? H , (13) 
with the inverse temperature (3 = 1/T, the partition function 

Z = exp(l + A ) = Tr e~ m , 

and the grand Hamiltonian 

H = H — p N - ^Ai - A . (14) 

Let us explain in more detail the role played by the Lagrange multipliers po and p\. 
Above the Bose-Einstein condensation point, when N Q = and Ni = N, one needs just one 
Lagrange multiplier p 1 = p, coinciding with the system chemical potential, whose role is to 
preserve the total average number of particles N =< N >. However, below the condensation 
point, when the global gauge symmetry is broken, there are two types of particles, condensed 
and uncondensed ones. The number of condensed particles Ao, according to the Bogolubov 
theory [8,9,17,18], has to be such that to make the system stable by minimizing the thermo- 
dynamic potential. For an equilibrium statistical system with the statistical operator (13), 
the grand thermodynamic potential is 

Q — —T In Tre -/W , (15) 



4 



with the grand Hamiltonian (14). Extremizing the grand potential (15) with respect to the 
number of condensed particles, from the equation 

m 



dN 
one obtains 

dH 

"° = < dN > ■ 

This means that the Lagrange multiplier po is responsible for the thermodynamic stability of 
the system. Another Lagrange multiplier, p±, guarantees the normalization condition (5) for 
the number of uncondensed particles Aq. But the latter, since Aq = N — N , implies that p 1 
preserves the total average number of particles N. In this way, for a Bose-condensed system, 
contrary to the system without the Bose-Einstein condensate, there are two conditions on 
the number of particles. One condition, as earlier, is that the total number of particles be 
N. And another condition is that the number of condensed particles, iVo, would be such 
that to provide the stability of the system by minimizing the thermodynamic potential. 
This is why one needs two Lagrange multipliers in order to guarantee the validity of these 
two conditions at each step of any calculational procedure. As is shown below by practical 
calculations, the use of two Lagrange multipliers makes the theory self-consistent, avoiding 
the Hohenberg-Martin dilemma and yielding the results that are in agreement with those 
derived analytically for the weak-coupling limit as well as obtained by Monte-Carlo computer 
simulations for strong interactions. 

It is also important to keep in mind that the introduction of Lagrange multipliers is 
a technical method allowing us to simplify calculations. The number of the introduced 
multipliers is connected with the concrete properties of the employed approach. Thus, in the 
Bogolubov theory [8,9,17,18] one deals with two independent field variables, the condensate 
wave function r)(r,t) and the operator of uncondensed particles ipi(r, t). This is why, as is 
explained above, it is convenient to introduce two Lagrange multipliers. 

One could ask whether we could limit ourselves by introducing a sole Lagrange multiplier. 
The answer is straightforward: Yes, we could, but the calculational procedure should then 
be changed. For instance, we could follow the way of Hugenholtz and Pines [22], which 
was also used by Gavoret and Nozieres [23]. They consider a uniform equilibrium system at 
zero temperature, defining the number of condensed particles N = N (p,T) as a function 
of density p and temperature T from the extremization of the internal energy E =< H >. 
The found N (p,T) is substituted explicitly into the Hamiltonian H, after which one works 
with the grand Hamiltonian H = H — pN±, where Aq is the operator for the number of 
uncondensed particles. Since the number of condensed atoms has already been defined earlier 
from the stability condition, one requires now to use the sole Lagrange multiplier aiming at 
guaranteeing the normalization condition Aq =< Aq > for uncondensed particles, hence, 
because of the fixed relation N\ = N — Nq, preserving the total number of particles N. This 
way of calculations, with details expounded in Refs. [22,23], is mathematically equivalent to 
the procedure, when No has not been fixed in advance but, instead, a Lagrange multiplier p Q 
is introduced to guarantee the stability condition in the process of calculations, after which 
the number of condensed particles is defined as Nq = N — Ni, which specifies the condensate 
depletion. It is this method of using an additional Lagrange multiplier that is employed in 
the present paper. 
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Even more, it could be possible to work introducing no any Lagrange multipliers at 
all, but which again would necessitate a different calculational procedure. Thus, one could 
use the Girardeau-Arnowitt approach [24,25] in the frame of the canonical ensemble, which 
requires to invoke the so-called number-conserving field operators. The weak point of this 
approach is that, as is well known [24,25], it yields the unphysical gap in the spectrum, 
when approximate calculations are involved. Girardeau mentioned [26] that an exact theory 
should not have the gap. Later, Takano [27] demonstrated that, really, the gap should not 
arise when all terms of the Hamiltonian are taken into account. Unfortunately, there are no 
exact solutions for a realistic interacting Bose system, so that one always has to resort to 
some approximations, in the course of which the gap again reappears. This is why it is more 
convenient to work with the grand canonical ensemble, introducing the Lagrange multipliers 
that would assure the self-consistency of any calculational scheme. 

Thus, the method of Lagrange multipliers has to be treated as a technical procedure 
allowing us to simplify calculations. The number of the required multipliers is intimately re- 
lated to the concrete calculational details. In the case of the Bogolubov approach [8,9,17,18], 
introducing two different field variables, the condensate function 77 and the field operator of 
uncondensed particles ifj±, it is convenient to define two Lagrange multipliers that guarantee 
the validity of two conditions, the thermodynamic stability condition and the conservation 
of the total average number of particles. 

When two Lagrange multipliers ate involved, none of them plays the role of the system 
chemical potential. To define the latter, we may proceed as follows. Keeping in mind that 
in experiments the total number of particles is usually fixed, we may write for the internal 
energy (11) the standard relation 

E = <H> + fiN , (16) 

connecting E with the average of the grand Hamiltonian < H > and with the system 
chemical potential /i. At the same time, substituting into Eq. (11) the grand Hamiltonian 
(14), and taking into account condition (8), we come to the expression 

E = < H > +fi N + fuNi . (17) 

Comparing Eqs. (16) and (17) gives the definition of the system chemical potential 

li = ii n + (18) 

expressed through the Lagrange multipliers /io and ji\ and the related atomic fractions 

The equations of motion for the variables 77 and tpi are given in the usual manner as 

■ d . , 5H 

i — ri(r,t) = -— 19 

dt /v ; Sr)*(r,t) K ' 

for the condensate function and, respectively, as 

at 5ipl(r,t) 
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for the field operator of uncondensed atoms. 

In this way, the representative statistical ensemble for an arbitrary equilibrium Bose sys- 
tem with broken gauge symmetry is the pair {^(ipi), p} of the Fock space of microstates 
and the statistical operator (13) with the grand Hamiltonian (14). The notion of rep- 
resentative statistical ensembles can be extended to nonequilibrium systems by considering 
the extremization of an effective action functional [28]. 



3 Bose Gas with Contact Interactions 

To specify the consideration, let us take the Hamiltonian energy operator in the usual form 
H = J ^t( r ) (- ^( r ) dv+ l -%J ^(t)^(t)^(t)^(t) dr ; ( 2i) 

corresponding to the contact interaction potential with the strength 

$ = 4tt - , (22) 
m 

where a s is the scattering length and m, mass. Here ^(r) = ^(r, t) is the shifted field 
operator (1). 

The evolution equation for the condensate function is obtained by averaging Eq. (19). 
To this end, we need the notation for the local condensate density 

Po (T,t) = \v(r,t)\ 2 , (23) 

normal density of uncondensed atoms 

Pl (r,t) =<4(r,t)i> 1 (*,t)> , (24) 

anomalous density 

<7!(r,f) =<Vi(r,0V>i(M)>, (25) 

and the triple correlator 

£(r,f) =<^(r,t)^ 1 (r,t)^i(r,t)> . (26) 

The total local density is 

p(T,t)=p (T,t)+ Pl (T,t) . (27) 

Then, averaging Eq. (19) yields the evolution equation for the condensate function 

+$ [p(r, t)r)(r, t) + Pl (r, t) V (r, t) + a ± (r, t)r)*(r, t) + £(r, t)} . (28) 

The equation of motion for the field operator of uncondensed atoms follows from Eq. 
(20) giving 



+$o [2p (r, t)^(r, t) + r?{r, t)4(r, t) + X(r, t) , 
where the last term is the correlation operator 

X(r, t) = [2^}(r, t)r)(r, t) + 77* (r, t)^(r, t) + 4(t, t)^(r, *)] ^i(r, t) . 

In equilibrium state, one has 



(29) 



0t 



77(r,t) = 0. 



(30) 



Also, if the system is uniform, then 



\q(r,t)\ 



No 
V 



= Po , 



piM) 



V 

N 



Pi , 



<ri(r,t) = (Ji , f(r,t) = f, p=_ = Po + Pl . 
In that case, equation (28) for the condensate function gives the Lagrange multiplier 



Po= [P + Pi + 0i + 



$ . 



(31) 



(32) 



The field operator of uncondensed atoms can be expanded over a complete basis. In 
general, if the system were nonuniform, in the presence of an external potential, it would be 
convenient to take a basis formed by natural orbitals [29,30]. For the uniform system under 
consideration, the natural orbitals are just plane waves v?fc( r ) = e' t]l ' r /\/V. The corresponding 
expansion of ip\ reads as 

^i(M) = J2 a k(t)<p k (r) . 

fc^O 

With this expansion, the grand Hamiltonian (14) takes the form of a sum 



H = J2 H (n) 



(33) 



?1=0 



of five terms, classified according to the number of the operators a k or a\ in the products. 
The zero-order term does not contain ak, 



(34) 



Generally, in order to satisfy condition (8), it is necessary and sufficient [31] that the 
Hamiltonian (14) would not contain the terms linear in ipi or a&. This can be achieved by 
choosing the corresponding Lagrange multipliers A(r, t) in Eq. (9). For a uniform system, 
because of the orthogonality condition (2), one has A = and = automatically. 

The second-order term is 



H ^ = £ 

k^O 

In the third-order term 



2m 



+ 2p $o -Aii J flfcO* + ^ Po$o (<4 a -k + a -kak) 



H {3) = $q (a\a q _ p a p + a\a\_ p a^ , 



(35) 



(36) 



p,q 
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the prime on the summation symbol implies that p ^ 0, q ^ 0, and p - q ^ 0. The 
fourth-order term is 

H(A) = ^7 2 Y! a i a \ a k+P a d-k , (37) 
ZK fc p,g 

where the prime means that p ^ 0, 0, k + p ^ 0, and k — q ^ 0. 

To proceed further, we need to invoke some approximation. A natural mean-field ap- 
proximation, in the presence of broken gauge symmetry, is the HFB approximation. This is 
used to simplify the higher-order terms (36) and (37). Employing the designations 

k 2 

LJ k = — + 2p$ - PL! (38) 

2m 

and 

A EE (p + <7i)$ , (39) 

we obtain 

_ - \ , , 

(40) 



H = E HFB + J2 ^fc a l°fc + it ( a l a -k + «-fc«fc) 



where the nonoperator term is 

E HFB = H^ - (2p? + a?) JV . (41) 

Hamiltonian (40) can be diagonalized by means of the Bogolubov canonical transforma- 
tion a k = u k b k + vl k b_ k , in which 

2 _ UJ k + E k 2 _ U k — £k 

Uk - ' Vk - ~^k~ ' 

with the Bogolubov spectrum 

e k = sfJl - . (42) 

Then one gets 

H = E B + Y,eAb k , (43) 

where _ 

E B ee E HFB + - J2(e k - cu k ) . (44) 

1 k^O 

By the Bogolubov [18] and Hugenholtz-Pines [22] theorems, the spectrum is to be gapless, 
which implies that 

lim e k = , e k > . (45) 

k— >0 

From here it follows that 

A*i = (p + Pi ~ <ti)$o ■ (46) 
The condensate multiplier (32) in the HFB approximation, when £ = 0, becomes 

Po = (p + Pi + ^i)$o • (47) 



9 



It is important to emphasize that the form of /xi in Eq. (46) is necessary and sufficient for 
making the spectrum gapless. That it is sufficient follows at once after substituting Eq. (46) 
into the Bogolubov spectrum (42). And the necessity stems from the Bogolubov theorem 
[18] which can be formulated as the inequality 

| /Jl --^ + E 12 (k,0)-E 11 (k,0)| < 7 '" 



2m ' 2m ' 

where S 12 and S n are the normal and anomalous self-energies. Setting here k = leads to 
the Hugenholtz-Pines relation 

// 1 = E 11 (0,0)-E 12 (0,0) . 

In the HFB approximation, one has 

S n (0, 0) = 2p% , £ 12 (0, 0) = (p + <7i)$ , 

from where one immediately obtains Eq. (46). 

With multiplier (46), spectrum (42) takes the form 



in which the sound velocity is 



k 2 V 



C EE J- . (49) 

V m 

For the diagonal Hamiltonian (43), it is straightforward to calculate all averages, such as 
the normal average 

n k ee < a\a k > (50) 

and the anomalous average 

°k = < ata-h > ■ (51) 
Their integration over momenta gives the density of uncondensed particles 



and, respectively, the anomalous average 



pl = J nk (W (52) 



dk 

0"! 



f dk , , 



The quantity \a±\ can be interpreted as the density of pair-correlated atoms [19]. 

In what follows, we concentrate on the zero-temperature properties of the system. When 
T = 0, then < b\b k >= 0. The normal average (50) is 

n k = — , (54) 
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while the anomalous average (51) becomes 

A . . 

°k = - 77— • (55) 

2S k 

Combining Eqs. (38), (39), (46), and (49), we have 

k 2 

uj k = h A , A = mc 2 . 

2m 

The equation for the sound velocity (49) can be represented as 

mc 2 = (p + <7i)$o ■ (56) 

Note that from here, in the limit of asymptotically weak interactions, we get the Bogolubov 
expression 

c^J^ (4>„-0). 
V m 

For the density of uncondensed atoms (52), we find 

fuj k -e k dk (mc) 3 

* = / -^r ^ = -3^- ■ (57) 

And for the anomalous average (53), we have 

A f 1 d k 

'' = -2/^(2^- < 58 > 

The integral / dk/Ek in Eq. (58) is ultraviolet divergent. To overcome this, we can re- 
sort to the standard procedure of analytic regularization [5,32]. For this purpose, we, first, 
consider the integral in the limit of asymptotically small <E>o, when the dimensional regular- 
ization is applicable, and then analytically continue the result to arbitrary interactions. The 
dimensional regularization gives 



/ 



1 dk 2 3/2 



e k (2tt)3 n 2 
Using this in Eq. (58), we obtain 



- — in- ~ \ />vA\) ■ 



0"! 



(mc) 2 

7l 2 



\Jmp $ . (59) 



It is important to stress that the anomalous density (59) enjoys the natural limiting 
property 

^i^O (po^O). (60) 

The physics of property (60) is evident. The existence of both the condensate density p 
and the anomalous density <j\ is caused by the gauge symmetry breaking. Both of them 
are nonzero as soon as the symmetry is broken, while both become zero if the symmetry 



11 



is restored. Any of these quantities could be treated as an order parameter for the broken- 
symmetry phase. So, both these quantities, p an d Ci, have to nullify simultaneously, when 
one of them tends to zero. 

We may also note that simplifying o\ by replacing ^mpo&o by mc, as is done in Ref. [14], 
is admissible only in the limit of weak interactions, when p Q ~ p, But for strong interactions, 
when p 0, this replacement does not hold, since then the limiting property (60) is not 
satisfied. Therefore, the results of Ref. [14] are quantitatively correct in the limit of weak 
interactions, though for strong interactions, they may give only a qualitative picture. And 
our aim in the present paper is to give a careful analysis of the system properties for arbitrary 
strong interactions in the whole range of $ £ [0, oo). This requires to employ Eq. (59), 
which explicitly satisfy the limiting condition (60). 



4 System Characteristics at Varying Interactions 

The effective interaction strength can be characterized by the dimensionless gas parameter 

7 = P 1/3 a s . (61) 

One often uses the quantity pa 3 as a parameter quantifying the interaction strength. How- 
ever, parameter (61), to our mind, is more convenient, since it better distinguishes between 
weak and strong interactions. Thus, the majority of experiments with ultracold trapped 
gases [1-7] deals with weakly interacting atoms, so that pa 3 s ~ 10~ 8 — 10~ 4 , which corre- 
sponds to 7 ~ 10 -3 — 10~ 2 . The effective interactions can be noticeably strengthened by 
loading atoms in optical lattices [33] . Contrary to weakly interacting gases, superfluid 4 He is 
a strongly interacting system. The 4 He atoms can be represented by hard spheres of diam- 
eter a s [34,35]. At saturated vapour pressure, one has [35] a s = 2.139A and p ks 0.022A~ 3 , 
hence pa 3 0.215, which is yet much less than one. But the corresponding 7 ps 0.599 is 
close to one. So, weak interactions are characterized by 7 1, while a strongly interacting 
system has 7 ~ 1. 

It is convenient to introduce the dimensionless sound velocity 



the fraction of uncondensed atoms 



pi/3 



ri! = — = — , (63) 



p N 
and the anomalous fraction 

a = ^ . (64) 
P 

In these dimensionless quantities, Eq. (56) writes as 

s 2 = 47r(n + a)7 . (65) 

Here 



n = 1 — n\ (66) 
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and 



while the anomalous fraction (64) is 

2s 2 

° = ^372 V^~o ■ (68) 

Four quantities, s, no, rii, and a, are defined by the system of four equations (65) to (68). 

In the limit of asymptotically weak interactions, when 7—^0 and o — > 0, the condensate 
fraction and sound velocity tend to the Bogolubov expressions 

n B = l-£= 7 3/2 (69) 

and, respectively, 

s B = ■ (70) 

For the higher orders with respect to 7, we find 

8 3/2 64 3 640 9/2 (n-\ \ 

« r- 1/2 16 2 32 7/2 3904 - 

^ 2 ^^ /+ TT + V5 T 2^ 7 ' (72) 

„ ^ * 7 3/2 + 32 3 _ M , /2 . (73) 

y 7T 7T 71-^ 

For asymptotically strong interactions, when 7 — > 00, we obtain 

n o-~ 7 3 - — hr 7 5 , (74) 



64 ' 512 V 9 J 

^(9^ _ » W^yV+^YV (76) 

4 [ 64 r 128 V 3 J 1 512 ^ 9 J 1 V ; 

So, the condensate fraction tends to zero, together with the anomalous fraction, as 7 — > 00, 
in agreement with condition (60). 

For the whole region of 7, we solve numerically the system of Eqs. (65) to (68), showing 
the results in Figs. 1,2, and 3. 

Figure 1 presents the condensate fraction (66) and this fraction (69) in the Bogolubov 
approximation. At small 7, up to 7 ~ 0.1, n and ub practically coincide with each other. For 
7 > 0.1, the Bogolubov approximation overestimates the condensate fraction. But % = 
at 7 = 0.762, while n is yet finite, though small. The condensate fraction of a homogeneous 
Bose gas, at zero temperature, as a function of the gas parameter was calculated by means 
of the Monte Carlo simulation by Giorgini et al. [36] up to 7 ~ 0.5. For this region of 7, 
our no is in good agreement with the Monte Carlo calculations. Monte Carlo techniques 
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have also been used for studying the condensate fraction of trapped atoms [37-39] . But the 
latter results cannot be directly compared with Uq in a homogeneous gas, since in traps, 
the condensate fraction is a function of spatial variables as well as of the trap shape [37- 
41]. What is possible and interesting to compare is the condensate fraction in superfluid 
4 He at zero temperature and our n at 7 ~ 0.6 corresponding to liquid 4 He with hard-core 
interactions. For 7 ps 0.6, we have Hq ~ 0.15, which is close to the condensate fraction 
no ~ 0.1, measured in experiments (as is reviewed in Refs. [3,42]) as well as obtained by 
Monte Carlo simulations (see review articles [43,44]). 

Figure 2 shows the dimensionless sound velocity (62) and its Bogolubov approximation 
(70). These quantities practically coincide up to 7 0.1. For 7 > 0.1 the Bogolubov form 
sb underestimates s till 7 ~ 0.7, after which it overestimates the latter. 

In Fig. 3, we compare the fraction of uncondensed atoms (67) with the anomalous fraction 
(68). As is seen, a > ri\ up to 7 ~ 0.7. The anomalous fraction a becomes substantially 
smaller than ri\ only for very large 7 ^> 1. This is in agreement with other calculations [45] 
confirming that anomalous averages cannot be neglected at low temperatures. 

Let us now analyse the ground-state energy 



The system chemical potential, defined in Eq. (18), is expressed through the Lagrange 
multipliers (46) and (47) and the fractions n and ni, which gives 



E =< H > +fiN 



(T = 0) . 



(77) 



According to Eq. (43), 



< H > = E B 



(T = 0) . 



H — (1 + ni + a — 2an 1 )p& . 



(78) 



For the dimensionless chemical potential 



p 2/3 



2mu 



(79) 



we get 



u = 8:^7(1 + rii + a — 1on\) . 



(80) 



From Eq. (44), we have 




(81) 



where the integral is calculated invoking the dimensional regularization [5], giving 




Equations (34) and (41) yield 




(82) 
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Summarizing these formulas, we find 



E p$ / o 9 \ 8(mc) 5 . . 



AT 2 V 1 / f57r 2 mp 

It is convenient to define the dimensionless ground-state energy 



E ° = -pmN' < 84 > 



for which we obtain 



E = 4tt 7 (l + n\ - 2an, - <x 2 + J^-) . (85) 
This can be compared with the Lee-Huang- Yang approximation [46-48] 

E LHY = 4tt 7 ^1 + 7 3 ^ , (86) 

derived for small 7. 

In the limits of asymptotically small and large 7, the dimensionless chemical potential 
(80) tends to 

95fi 

p _, 8 7T7 + — v^F 7 5 / 2 (7 -> 0) (87) 
3 

and, respectively, 

77 — * 16717 (7 — > 00) . (88) 
The dimensionless ground-state energy (85) at small 7 1 possesses the expansion 

E ~ 47T 7 + ^- v^F 7 5 / 2 + ^- 7 4 , (89) 
to y 

which reproduces the Lee-Huang- Yang approximation (86) for 7 — > 0. And for large 7 ^> 1, 
we find 

£„ ^8. 7 + 5 (9^) 1/3 - * (3^) 1/3 T-' + l (3,") 1/3 7- 4 . (90) 

Note that expressions (87) and (88) can also be obtained from Eq. (89) and (90) using the 
relation /i = dE/dN, valid for T = 0. 

Figure 4 illustrates the behaviour of the ground-state energy (85) and the Lee-Huang- 
Yang approximation (86). The latter is known to practically coincide with the energy calcu- 
lated trough the Monte Carlo simulations [36] up to 7 0.4. As is seen from the figure, our 
E is also very close to E LHY in the region < 7 < 0.4, but is lower than E LHY for 7 > 0.4. 
Hence, E well reproduces the available data of Monte Carlo calculations up to 7 ps 0.4. 



5 Conclusion 

The notion of representative statistical ensembles is applied to Bose systems with bro- 
ken global gauge symmetry. A general procedure is described for constructing the grand 
Hamiltonian for the representative ensemble of an arbitrary equilibrium Bose system. A 
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self-consistent mean-field theory is developed, which is both conserving and gapless. The 
properties of a uniform Bose gas at zero temperature are studied both analytically and nu- 
merically for the gas parameter varying between zero and infinity. Thus, in the frame of 
the suggested approach, strongly interacting systems can also be considered. For instance, 
as is known [35,43], some of the properties of superfluid 4 He can be understood by treating 
the potential as a hard-core interaction of diameter a s = 2.139A, which, at saturated vapor 
pressure, corresponds to 7 ~ 0.6. For the latter 7, we find the condensate fraction no of order 
10%, which agrees with the condensate fraction in helium at zero temperature, measured 
in experiments [42] and found in Monte Carlo simulations [43], being also of order 10%. 
Another application of the developed self-consistent mean-field theory with arbitrary strong 
interactions could be the description of Bose-Einstein condensation of multiquark clusters 
in nuclear matter [49,50]. 
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Figure Captions 



Fig. 1. Condensate fraction no (solid line) and its Bogolubov approximation % (dashed 
line) as functions of the gas parameter 7. 

Fig. 2. Dimensionless sound velocity s (solid line) and its Bogolubov approximation sb 
(dashed line) as functions of the gas parameter 7. 

Fig. 3. Fraction of uncondensed atoms n\ (dashed line) and anomalous fraction a (solid 
line) as functions of the gas parameter 7. 

Fig. 4. Dimensionless ground-state energy E (solid line) and the Lee-Huang- Yang 
approximation E LHY (dashed line) as functions of the gas parameter 7. 
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Figure 2: Dimensionless sound velocity s (solid line) and its Bogolubov approximation 
(dashed line) as functions of the gas parameter 7. 



21 




Figure 3: Fraction of uncondensed atoms n\ (dashed line) and anomalous fraction a (solid 
line) as functions of the gas parameter 7. 
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Figure 4: Dimensionless ground-state energy E (solid line) and the Lee-Huang- Yang 
proximation E LHY (dashed line) as functions of the gas parameter 7. 
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